Off-equilibrium fluctuation-dissipation relations 
in the 3d Ising Spin Glass in a magnetic field. 
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We study the fluctuation-dissipation relations for a three dimensional Ising spin glass in a magnetic 
fleld both in the high temperature phase as well as in the low temperature one. In the region of 
times simulated we have found that our results support a picture of the low temperature phase with 
broken replica symmetry, but a droplet behavior can not be completely excluded. 
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I. INTRODUCTION 

The understanding of the behavior of a spin glass in 
a magnetic field is a challenging issue from both experi- 
mental and theoretical sides. 

In the theoretical side, there are two competing theo- 
ries. In the droplet model the spin glass phase is unsta- 
ble (for any amount of magnetic field) and so there is no 
phase transition: there is only a pure state describing all 
the Ising spin glass in a magnetic field On the other 
side. Mean Field (MF) predicts a phase transition be- 
tween two phases. The first one is characterized by one 
pure state, and the low temperature phase is described 
by a countable number of pure states. In the Mean Field 
approximation a third order phase transition^ has been 
found between those two phases^ separated by the de 
Almedia-Thouless line See references for a descrip- 
tion and0 for a critic of the RSB picture. 

Hence, those two competing theories have opposite 
predictions about the overall behavior of a spin glass in 
the presence of a magnetic field. However, to perform 
experimental or numerical tests of the previous analyti- 
cal predictions has proved very diSicult despite the clear 
theoretical predictions (phase transition or not!). 

A further step in Mean Field computations is to take 
into account the effect of fluctuations. This can be done, 
for instance, using Field Theoretic methods and has been 
done in the past. Working with a projected theory 
(by taking only the replicon sector which contains the 
most divergent terms of the initial Hamiltonian) no fixed 
points have been found in the model. ^ It is important to 
mention that the existence or not of a transition in mag- 
netic field affects the existence or not of a phase with 
replica symmetry breaking (RSB) at zero magnetic field. 
In particular the absence of a fixed point for the Ising 
spin glass in a magnetic field supports the droplet pic- 
ture against the Mean Field one. 

A recent and detailed analysis^ reaches the same con- 
clusions of Ref. 0, and so three possibilities are opened: 



i) no phase transition at all, ii) first order phase transi- 
tion driven by fluctuations and finally Hi) second order 
phase transition dominated by a fixed point outside the 
accessible perturbative region (i.e. region of small param- 
eters). However, a fourth possibility has been recently 
opened by Temesvari and De DominicisjS who have ex- 
tended the field theoretic analysis further. They analyze 
a field theory in which the replicon and anomalous sec- 
tors are both critical going beyond the old analysis where 
only the replicon sector was taken critical. The main re- 
sult of this analysis is the existence of a new critical point 
(which appears at eight dimensions) taking the control of 
the phase transition at six dimensions. The authors also 
pointed out that this new fixed point provides a phase 
transition which has different features to that of Mean 
Field. 

On the numerical side the situation is a bit clearer (but 
not enough!). Looking at the off equilibrium numerical 
simulations, the difference between the mean overlap and 
the minimum overlap has been computed^Siiiii^ In four 
dimensions and not for too low temperatures there is a 
clear difference between these two measurements, which 
is a clear signature for Replica Symmetry Breaking. An- 
other off equilibrium approach is to compute the violation 
of the fluctuation dissipation theorem out of equilibrium. 
This has been done using a slightly modified version of 
the four dimensional Gaussian spin glass, and it was 
found that the violation is well understood in terms of a 
non trivial low temperature phase. In this paper we will 
follow this approach but working in three dimensions and 
simulating the Edwards- Anderson model. It is important 
to notice, that the same kind of studies on the violation 
of the fluctuation-dissipation out of equilibrium can be 
done in real experiments. In fact, in a recent paperji^ 
the violation of the fluctuation-dissipation relations in 
an Ising spin glass (in zero magnetic field) has been re- 
ported, and the experiment can be explained in terms of 
Replica Symmetry Breaking. 

Another numerical method is to use exact ground state 
techniques in order to understand the qualitative features 
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of the low temperature phase. This has been done in 
three dimensions and a RSB behavior has been found 
between zero and a magnetic threshold (for the Gaus- 
sian Ising spin glass in 3d, this threshold is near 0.65), 
however this work cannot exclude completely a droplet 
behaviorpi^ A recent studj»i^ points out that there is no 
phase transition at all (at finite temperature) but they 
cannot exclude a critical field below 0.4. 

We can cite different numerical studies working at equi- 
librium. Those studies have been mainly done in four di- 
mensions in order to avoid the proximity to the lower crit- 
ical dimension (which is between two and three, at least 
in the Ising spin glass without magnetic field). Those 
studies are not fully conclusive, but the existence of a 
finite temperature phase transition emerge as the most 
likely explanation of the numerical data.^^'^^ 

On the experimental side the situation is not clear. 
There is strong experimental evidence about an irre- 
versibility line (where the zero field cooled (ZFC) and 
the field cooled (FC) magnetization start to be different), 
but unfortunately, that line depends on the time, and 
so, we have only off equilibrium information about what 
happens in the presence of a magnetic field. In the liter- 
ature, one can find some attempts to analyze the scaling 
behavior of the freezing temperatures and the conclusion 
is that no phase transition exists»i£ 

However, recent experimental studies based in the 
fiuctuation-dissipation relations point out the existence 
of a phase transition in the presence of magnetic fieldii^ 
Moreover, a phase transition has been reported in a 
Heisenberg spin glass (AuFe) in three dimensions in the 
presence of a magnetic field against the droplet predic- 
tion (we recall that the droplet model predicts no phase 
transition independently of the number of components of 
the spin)>S& 

We will study the three dimensional Ising spin glass 
using an off-equilibrium approach based on the computa- 
tion of the fiuctuation-dissipation relations. This method 
has provided and important tool to investigate the low 
temperature properties of disordered systems (and it has 
been very useful in the study of non disordered systems 
such as glasses) iSiiSSiSiSi 
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By < ij > we denote the sum over nearest neighbor pairs. 
The Jij are chosen from {-1-1, —1} randomly and h is the 
external magnetic field. We have studied systems with 
magnetic field h = 0.2 and lattice sizes L = 20, L = 30 
and L — 60 

We have used the SUE parallel computer. SUE is a 
dedicated machine with an overall performance of 0.22 
ns per spin flip. See references i25u26i for a detailed de- 
scription of this computer. 

Given a quantity A{t) that depends on the local vari- 
ables of our original Hamiltonian (7i), we can define the 
following autocorrelation function 



C{h,t2) = {A{h)A{t2)) 
and the response function 

S{A{h)) 



R{ti,t2 



SAh{t2) 



(2) 



(3) 



Ah=0 



where we have assumed that the original Hamiltonian 
has been perturbed to 



Ah{t)A{t)dt . 



(4) 



The brackets (• • •) in eq. |5J and eq. (0) imply here a 
double average, one over the dynamical process and one 
over the disorder. 

As usual one could choose A{t) = ai{t) and the re- 
sponse function should be 



R{ti,t2 



S'm{ti) 



SAh{t2) 



Ah=0 



6Ah{t2) 



(5) 



Ah=0 



where m{t) = (o'i(O) 

However, to improve the signal of the autocorrelation 
we have used in the present paper: 



II. THEORETICAL BASIS 

We have focused this paper on the study of the fluctu- 
ation dissipation relations in the off-equilibrium regime. 
To do this we need to define the spin-spin autocorrela- 
tion function and the response of the magnetization to a 
small change of the magnetic field of the system. 

In order to make the paper self-contained and to fix the 
notation, we shall recall some important results about the 
off equilibrium fluctuation-dissipation relations. We have 
simulated the binary Ising spin glass in three dimensions 
on a cubic lattice of volume V = L'^ with helical bound- 
ary conditions. The Hamiltonian of the system is given 



1 y 

C{tiM) = -^(fT,(ii)a,(t2)) , 

i—1 



(6) 



and m{t) ~ jF^iWiit)) ■ We remark that we are inter- 
ested in global fluctuation dissipation relations. Recently 
work has been done in local microscopic fluctuation dis- 
sipation relations)^ but we will not study them in this 
paper. 

In the dynamical framework, assuming time transla- 
tional invariance, it is possible to derive the fluctuation- 
dissipation theorem (FDT), that reads 



R{hM) = (iOih - t2) 



dt2 



(7) 
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where the inverse temperature is /3 = 1/T. 

The fluctuation-dissipation theorem holds in the equi- 
librium regime, but in the early times of the dynamics we 
expect a breakdown of its validity. Mean Field studies^S 
suggest the following modification of the FDT: 



R{h,t2) ^ f3X{ti,t2)e{h -t2) 



dC{h,t2) 
dt2 



(8) 



It has also been suggested inSSiSSiffi that the func- 
tion X{t,t') is only a function of the autocorrelation: 
X{t,t') = X{C{t,t')). We can then write the follow- 
ing generalization of FDT, which should hold in early 
times of the dynamics, the off-equilibrium fluctuation- 
dissipation relation (OFDR), that reads 



R{h,h)=l3X(C(t,,h})Hti 



Ot2 



We can use the previous formula, eq. (O, to relate the 
observable quantities defined in eq. l(2Jl and eq. Using 
the functional Taylor expansion we can write 



m[ft.-f A/i](t) = m[h]{t) 



' 8m[h']{t) 
Sh'{t') 



Ah{t') + 0{Ah^) , 



(10) 



h'{t)=h{t) 



and so. 



Am[/i, A/i](t) = I dt' R{t,t')Ah{t') + 0{Ah^) . (11) 

where we have defined Ani[h, Ah]{t) = m[h + Ah]{t) — 
m[h]{t). Eq. (|ll|l is just the linear-response theorem 
neglecting higher orders in Ah. By applying the OFDR 
we obtain the dependence of the magnetization with time 
in a generic time-dependent magnetic field (with a small 
strength), Ah{t), 



Am[/i,A/i](t) ~/3 



* dt' X[Cit,t')]^^^AHt') 



(12) 

Next we let the system evolves with the unperturbed 
Hamiltonian of eq. ^ from t = to t = t^, the so 
called waiting time, and then wc turn on the perturbing 
magnetic field Ah (hence, the system feels a magnetic 
field h + Ah). Finally, with this choice of the magnetic 
field, we can write (ignoring in our notation the fact that 
Am also depends on t„) 

Am[h,Ah]{t) ~ A/i/3^ dt' X[C{t,t') f'^^^P , (13) 

and by performing the change of variables u = C{t,t'), 
eq. ((T^ reads 



Am[h, Ah]{t) ~ Ahf3 [ du X[u] 



(14) 



where we have used the fact that C{t, t) = 1 (always true 
for Ising spins). In the equilibrium regime (X = 1 as the 
fluctuation-dissipation theorem holds) we must obtain 



Am[h,Ah]{t) ~ A/i/3(l-C(i,tw)) 



(15) 



i.e. Am[h, Ah]{t)T/ Ah is a linear function of C{t,t^) 
with slope —1. We remark that we can use this formula 
to obtain qmax as the point where the curve Am[/i, Ah]{t) 
against C{t,t^) leaves the line with slope —13 Ah. 

In the limit t,t^ oo with C{t, t^) — q, one has that 
X{C) x{q), where x{q) is given by 



x{q) 



dq' P{q') 



(16) 



where P{q) is the equilibrium probability distribution of 
the overlap. Obviously x{q) is equal to 1 for all q > 
9max, and we recover FDT for Cit^t^) > gmax- This 
link between the dynamical function X{C) and the static 
one x{q) has been already verifled for flnite dimensional 
spin glasses i2S The link has been analytically proved for 
systems with the property of stochastic stabilityi^i 
For future convenience, we deflne 



S{C) = / dq x{q) 



or equivalently 



d^S{C) 



d^C 



In the limit where X 



C=q 

X we can write eq. H14|l as 



Am[Ah]{t) T 
'Ml 



S{C{t,t^)) 



(17) 



(18) 



(19) 



Looking at the relation between the correlation func- 
tion and the integrated response function for large we 
can thus obtain gmaxi the maximum overlap with non- 
zero P{q), as the point where the function S{C) becomes 
different from the function 1 — C. 

From the function S(C) we can get information on the 
overlap distribution function P{q), through eq. ((TH|) . Let 
us recall which is the prediction for the S{C) assuming 
the validity of each one of the competing theories de- 
scribed in the introduction. The droplet model predicts 
P{q) — 5{q — q) and consequently 



5(C) = 



- q for C < g , 
C for C> g . 



(20) 



In models with only one state, as the droplet model 
predicts for this model, the equilibrium time is flnite ir- 
respective of the value of the volume of the system, hence, 
we can always thermalize any volume, and so the asymp- 
totic behavior, for waiting times larger than the equili- 
bration time, is composed only for the straight line 1 — C. 
There is no horizontal part. 
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On the other hand the MF Hke prediction for the over- 
lap distribution^ P{<}) = (1 ~ XM)S{q — 9max) + ^M^iq — 
Qmin) +p{q) (where the support oi p{q) belongs to the in- 
terval [(jmin,gmax], gmin « h'^^^ and Qmax mainly depends 
on the temperature), implies that 



5(C) 




for C < Qmin , 
for gmin <C < Qn 

for C> 



(21) 



where s{C) is a quite smooth and monotonically decreas- 
ing function such that 



(f~s{C) 



dC2 



(22) 



C=q 



To finish this section we will recall an approximate scal- 
ing property of the probability distribution of the overlap 
that was introduced by Parisi and Thouless (hereafter 
PaT).'^^ In particular in Mean field the PaT hypothesis 
implies 



S{C) 



1 - C for C > 9,nax , 

TVT~C for q^in <C<qn 



(23) 



The result for C > gmax is general (and true for finite 
dimension) and for qmin < C* < ^max we make the follow- 
ing Ansatz: S{C) = AT{1 - C)^ (in Mean Field A = 1 
and B = 1/2). If we substitute this Ansatz in eq. IjlOfl 
we obtain the following scaling equation 



mT 



T-'t' = f {{1 ~ C)T-'t>) 



(24) 



where / is a scaling function and (j) = 1/(1 — 5) (in 
Mean Field (j) = 2). In order to be consistent the scaling 
function should be composed by a linear part (x) and by 
a power law part (Ax^). 

We have only measured the autocorrelation function 
(see eq. ©) and the response function (see eq. jSJ). 



III. NUMERICAL RESULTS 
A. On the critical temperature 

Assuming the existence of a phase transition in mag- 
netic field, we can estimate the shift of the critical tem- 
perature when a small magnetic field is turned on using 
the Mean Field approximation. The main formula isS 



Tc{0) 



1/3 



2/3 



(25) 



where J is defined (in Mean Field) as Jfj — ,P/N, N be- 
ing the volume of the system (or the coordination number 
in the Mean Field approximation) , being the random 



couplings between the spins, Tc(0) being the critical tem- 
perature in absence of magnetic field and finally Tc{h) be- 
ing the critical temperature in the presence of a magnetic 
field h. That is the formula that fixes the AT instability 
in infinite dimensional spin glasses. 

We can modify that formula (eq. (|25|l ) for a finite 
coordination number. Let z be the coordination number 
of our lattice. We recall that our have unit variance 
and so J = y^. Hence we can write 



Tc(O) 



1/3 



2/3 



(26) 



In our case z = 6, Tc(0) ~ 1.14, and so T{h = 0.2) ~ 
0.945 {h = 0.2 is the magnetic field simulated in the 
present work). Notice that near zero magnetic field 
the phase transition line has vertical slope {dT{h)/dh ~ 

l//ll/3). 

In order to check the existence or not of a phase tran- 
sition (using the OFDR as a tool) we have simulated at 
very high temperature (T — 2.5) and a lower tempera- 
ture (which is below our previous estimate of the critical 
one, T = 0.714). 



B. OFDR in the high temperature region 

We have simulated the system at temperature T — 2.5 
in a magnetic field h = 0.2 and with perturbing fields 
Ah = 0.01 and Ah = 0.03 (in order to check linear re- 
sponse) and different waiting times: 409600 and 819200. 
The number of samples simulated was about 6400 sam- 
ples for each waiting time. 

We show in figure ^ the plot for ty, ~ 819200 and 
L = 30. 

For the largest value of the waiting time simulated all 
the data stay on the equilibrium line 1 — C (i.e. this 
waiting time is greater than the equilibration time for 
this lattice size). 

In the paramagnetic phase, droplet and RSB agree: 
for a "finite" volume and very large times (grater than 
the equilibration time) all the points should lie on the 
straight fine Am{t)T / Ah = 1-C with C G [0, ^ea] (equi- 
librium situation). For intermediate situations (i.e. not 
so large waiting times) the curves lie below the straight 
line (see figure |2Jl and the final straight line, is built 
from below (i.e. curves with lower waiting times lie be- 
low those with higher ones). This behavior is similar 
to that found in the two dimensional Ising spin glass in 
a magnetic field for a finite temperature (the system is 
paramagnetic).'^^ See also figure 5 of reference^*' for an 
example of a FDT plot in a paramagnetic phase in the 
four dimensional Ising spin glass. 



C. OFDR in the low temperature region 

The situation at a lower temperature is dramatically 
different. 
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C(t,t ) 



FIG. 1: Off equilibrium fluctuation-dissipation relations for 
T = 2.5, h — 0.2 and L = 30. We have drawn the equilibrium 
straight line 1 — C. We plot AmT/h against C{t,t„) for the 
waiting time = 819200. 




FIG. 2: As figurelUbut with t„ = 409600. 



We start showing the numerical results for one of the 
lowest temperature simulated, T — 0.714. All the sim- 
ulations reported in this subsection have been done at 
/i = 0.2 as the external magnetic field. 

In figure|21we show AmT/A/i against C{t,t^) for dif- 
ferent waiting times iw and perturbing magnetic field 
Ah = 0.03 for the L = 30 lattice. 

The first check we have performed is to control that 
we are in the linear response regime. To do this we have 
computed the OFDR for different perturbing magnetic 
fields Ah = 0.01 and 0.03. We have found that the re- 
sults are independent of these two values of Ah. In the 
following we show the results obtained with Ah = 0.03. 

The second check has been to verify that our results 
are lattice size independent. In figure 01 we can see the 
results for L = 30 and L — 60, with perturbing field 
Ah = 0.03. For L = 20, the perturbing field Ah = 0.03 
has proved too noisy, and we have used Ah = 0.06. It 



can be seen that the behavior of i = 30 is asymptotic in 
this kind of simulations (i.e. for the time scales that we 
have simulated): the L — 20 points are still a bit noisy 
but L = 30 and L = 60 coincide. Using this information 
we will focus on the L = 30 lattice in the rest of the 
paper. We can state that we have simulated 512 samples 
for tw = 81920, 416 samples for = 163840 and 3232 
for t„ = 327680 and 1638400 in the L = 30 lattice. In 
addition 6200 samples in the i = 20 lattice and 412 in 
the largest lattice that we have simulated {L = 60). 

The third goal is to study the dependence on the 
waiting times of the OFDR curves. We have simulated 
t„ = 81920, 163840, 327680 and 1638400 and from figure 
13 it is possible to see that the curves rise as the waiting 
time is larger. Moreover, the curves for the larger waiting 
times are just compatible within our error bars. In this 
sense we are confident that the curve corresponding to 
tw — 1638400 represent very well the overall behavior of 
the system (very large volumes and times or equivalently 
infinite volume and waiting times) . This behavior is very 
important because our final curve (i.e. tw = 1638400) has 
a clear curvature which should be absent if the droplet 
model holds (see eqs. (|^ for droplet and (|^ for RSB 
predictions). 

We have plotted in figure O two additional straight 
lines. The first line, horizontal, corresponds to the 
asymptotic value of Am{t)T / Ah. To obtain this value, 
we have performed a simulation reaching times much 
longer than the ones used in the OFDR curves, but using 
a smaller number of samples. In figure [S] we can see the 
evolution of m{t). The horizontal lines are the error band 
related to the asymptotic value of Am{t) = 0.1286(15). 
For this value we do not need to do any kind of extrapo- 
lation since we have reached the thermodynamic value of 
Am{t)T/Ah in this long simulation (O(IO^) Monte Carlo 
steps). 

If the droplet model holds, the final (asymptotic) curve 
should be composed by the 1 — C straight line (C S 
[^EA,!]) as explained above. This implies that we can 
estimate the "droplet" prediction for the order parameter 
as qea — 0.694(4). We have marked this value with a 
vertical line in figure |31 

We would remark at the end of this section the follow- 
ing points: 

• We have obtained a tw-independcnt final curve (and 
L-independent) , at least within our statistical pre- 
cision. We believe that this curve represents with 
high accuracy the behavior for large volumes and 
times of an Ising spin-glass at T = 0.714 and 
h = 0.2. In this scenario we can estimate that 
qea — 0.76(2) (the points in which the points leave 
the straight line) which differs from the droplet 
value 0.694(4). 

• We cannot avoid a dependence on the waiting times 
beyond our numerical precision, and so we can not 
exclude completely a droplet phase with (/ea = 1 ^ 
Am{oo)T/Ah ~ 0.694(4). 
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C(t,t,) 



FIG. 3: Off equilibrium fluctuation-dissipation relations for 
T = 0.714, L = 30 and /i = 0.2. We have marked the equilib- 
rium straight line 1 — C. We plot AmT /h against C(t, t„) for 
three different waiting times. We have also plotted the error 
band for the asymptotic value of /S.mT /h. 
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FIG. 4: Off equilibrium fluctuation-dissipation relations for 
T — 0.714 and three different lattice sizes L — 60, L — 30, 
L = 20. The perturbing field is 0.03 for L = 30, 60 lattices 
and 0.06 for L = 20 lattice. We show the data for one of the 
larger waiting time simulated t„ — 327680. 

• As we have cited in the section devoted to high tem- 
perature, the final curve is built from below. At low 
temperature, in the droplet scenario, we should ex- 
pect the same behavior as at high temperature and 
so the final curve should build from below. The 
point is whether the curves for large waiting times 
stop or not before they reach the droplet prediction. 
Our numerical data suggest that the curves stop 
before the droplet final curve, and that the asymp- 
totic curve shows the characteristic curvature of a 
phase with RSB. 

We end this section showing a figure corresponding to 
the crossover region. In the following discussion we will 



■ llfc-jjl.ll J.lW.J ,to-iii.j-j^jTyjifai,^^fliji,jji] 
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E 
< 

0.005 - 
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le+07 2e+07 3e+07 4e+07 5e+07 6e+07 7e+07 

t 

FIG. 5: Am{t) against t for U = 327680 and L = 30 in a 
simulation longer in time than the ones used in the figures of 
FDT, but with much less samples 1152. The fit to equilibrium 
plateau is also shown (we show again the error band for the 
mean value). 

restrict ourselves to a qualitative level. In figure |3 we 
have shown the OFDR for the three values of tempera- 
ture (T = 1.25, 1.11 and 1.0) and t^, = 327680. It is clear 
that the largest temperatures shows a clear signature of 
a paramagnetic phase (i.e. small curvature and almost a 
straight line) . You can compare these curves with a clear 
paramagnetic one, see figure |21 which also shows a small 
curvature at the end of the curve. We remark that the 
critical temperature of the model with no magnetic field 
is about 1.14, and so we are (qualitatively) exploring the 
region near the vertical of this point. On the basis of 
mean field we should expect that the line of transitions 
emerges from the point at zero magnetic field with verti- 
cal slope. Hence, this plot and a transition temperature 
near 1.1 is compatible with this scenario. The prediction 
using eq. (|26|l was 0.95. Obviously, in order to see a clear 
paramagnetic curve with no violations, we refer to figure 
n We remark that this methods based in the violation 
of fluctuation-dissipation is not so powerful to determine 
with precision transition points; is a good method to de- 
cide if one point (T, h) (well inside the phase in order to 
avoid the crossover region) behaves in a way or not. 



D. OFDR in the high magnetic field region 

In this section we study the properties of OFDR at 
a fixed low temperature when the magnetic field grows. 
The goal of this section is to find when the behavior of the 
OFDR relations change from a non trivial one (as found 
for T = 0.714 and h = 0.2) to a trivial one (droplet) 
as far the magnetic field becomes larger. Notice that 
this temperature (T = 0.714) is far away of the critical 
temperaure of the model with no magnetic field {T — 
1.138), avoiding so, crossover effects between the phase 
transition at zero field. 
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C(t,t^) 

FIG. 6: We show OFDR for three temperatures, for t„ — 
327680 and L = 30. 



We show in figures [7| and |S1 the results obtained at 
T = 0.714 and h = 0.4 and h — 0.6 respectively. 

We start discussing the h — OA plot (figure [7)). If we 
compute, as above, gEA as the minimum value of the 
correlation (starting from C = 1) for which the points 
do not lie (using one standard deviation as critera) on 
the straight line, we obtain that this value of the field is 
still not statistically compatible with the droplet value. 
If we relax the one standard deviation criteria to two or 
three standard deviation, the behavior can be described 
as droplet. Obviously for larger magnetic field we can 
show analytically that the behavior is droplet (the mag- 
netic contribution in the Hamiltonian becomes dominant 
and we can drop the spin glass term). 

For a larger magnetic field, h — 0.6, the situation is 
clearer (figure |HJl. Practically all the points are on the 
straight line (slightly below but always at a distance less 
than a one- two standard deviations). We have a hori- 
zontal part (the system has reached its asymptotic mag- 
netization) since our largest waiting time is still smaller 
than the equilibration time for this magnetic field (in the 
droplet, this time is finite). 

The conclusion of this section is that h = 0.4 the situa- 
tion is still not clear but ft, = 0.6 is droplet. We have ob- 
served a clear change in the behavior of OFDR in the re- 
gion h ~ 0.4 — 0.6. Assuming a phase transition (between 
RSB and droplet) this impHes that hc{T = 0.714) ~ 0.6. 

This figure compares very well with that obtained at 
zero temperature in a related model (which uses Gaussian 
coupling). It has been obtained a critical magnetic field 
at zero temperature of 0.6i^ or 0.4 if we use reference 

m. 




C(t,t,) 

FIG. 7: Off equilibrium fluctuation-dissipation relations for 
T = 0.714, L = 30 and h = 0.4. We have marked the equflib- 
rium straight line 1 — C. We plot AmT/h against C{t, t„) for 
three different waiting times. We have also plotted the error 
band for the asymptotic value of AmT/h. 




C(t,t„) 

FIG. 8: Ofi^ equilibrium fluctuation-dissipation relations for 
T = 0.714, L = 30 and /i = 0.6. We have marked the equilib- 
rium straight line 1 — C. We plot AmT/h against C{t, t^) for 
two different waiting times. We have also plotted the error 
band for the asymptotic value of AmT/h. 



E. Scaling properties of OFDR in the low 
temperature region 

In this section we will study the scaling properties of 
the iw and L-independent fluctuation dissipation curves 
obtained for different temperatures. The main goal of 
this section is to study the degree of accuracy of the 
approximate PaT Ansatz in the three dimensional Ising 
spin glass in a magnetic field. That Ansatz has been 
found very adequate to describe the low temperature 
fluctuation-dissipation curves both in three and four di- 
mensions in the presence of a magnetic field. Moreover 
this scaling Ansatz has been checked in experiments find- 
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C(t.t^) 

FIG. 9: Off equilibrium fluctuation-dissipation relations for 
five different temperatures. We have marked the equilibrium 
straight lines. Notice that we plot Am//i instead AmT/Zi as 
in previous plots, hence the slopes of the equilibrium lines 
are l/T. We show only the data computed with the largest 
waiting time tw = 327680. 




C(t,t ) 



FIG. 10: Off equilibrium fluctuation-dissipation relations for 
T = 0.625 (the coldest temperature we have simulated) for 
i = 30, /i = 0.2 and A/i = 0.03. We show the data for 
the largest waiting time simulated t„ = 327680 and for — 
162840. The date seem to be asymptotic in the statistical 
error. 



ing that it describes very well the experimental data. The 
Ansatz has also been studied in the two dimensional spin 
glass (no phase transition) and it has been found that the 
curves computed for different temperatures at the same 
{t^ = 10"*) for a large lattice {V = 400^ ) also follow 
this Ansatz (see figure 1-a of reference . Notice that 
this waiting time is not an asymptotic one (see fig 1-b of 
reference l33^ . 

We show in figureElthe and L-independent fluctua- 
tion dissipation curves obtained at five different tempera- 
tures. In particular, we show in figure lTUl the fluctuation- 




(1/T)''(l-C(t.tj) 

FIG. 11: Scaling plot of the off equilibrium fluctuation- 
dissipation relations for four different temperatures in the low 
temperature region. 



dissipation plot for the lowest temperature we have sim- 
ulated, T — 0.625, and for two different and large wait- 
ing times (iw = 162840 and = 327680). It is clear 
that these two waiting times are asymptotic (within er- 
ror bars) and so we are again confident in that our curves 
plotted in figure El are asymptotic. 

We have deleted the temperature factor in the ordinate 
axis and so each FDT straight line has slope l/T. We 
can see that when the data leave the pseudo-equilibrium 
region (i.e. the straight line), they go to the same curve 
(independently of their temperatures). This is a strong 
signature that the PaT Ansatz works. 

We have tried in figure ^] this kind of Ansatz and 
the scaling is very good and so we are confident that the 
PaT scaling describes with great accuracy the behavior of 
the fluctuation dissipation curves in a large temperature 
window. 

Notice that the PaT scaling works for our L— and 
t„— independent curves. We have found a good scaling 
for values of S (1.2,1.4). In figure ITTI we use central 
value for (j) — 1.3. Two clear and distinctive regimes 
can be seen in that figure . The first one correspond to 
the quasi-equilibrium regime: in that part of the figure 
the behavior is linear and so it matches with the quasi- 
equilibrium regime AmT/h = 1 — C The second one 
corresponds to the aging regime: that part of the plot 
can be parametrized with a power law with the B expo- 
nent introduced above in the paper. We have obtained 
B — 0.27(3) which provides cj) — 1.37(6), which is a com- 
patible value with the (j) value used in the scaling plot 
(this is a check of consistency of the scaling law!). For 
completeness we report the value of A (we remark that 
the aging region follows a law Ax^ , where x is the scaling 
variable T-'I>{1 - C)): A = 0.52(1). 

We can compare the values obtained for A and B with 
previous results published in the literature. In the 3d 
Ising spin glass with no magnetic field A = 0.7 and B = 



9 



OAlM For the Ad Ising spin glass again with no magnetic 
field A = 0.52 and B = 0.41^2 

We can see that in absence of magnetic field the B 
value is close to the Mean Field value (0.5) whereas the 
magnetic field value in three dimensions is clearly far 
from the MF value. 

Following reference this kind of scaling it is not 
enough to detect a RSB phase (they found in the two 
dimensional Ising model — with no phase transition at fi- 
nite temperature — a PaT scaling for their OFDR) . Nev- 
ertheless, inS^ the PaT scaling only works for points with 
the same waiting time, instead, in our plot we have points 
computed with different waiting times. In effect, we re- 
mark again, our scaling is tw-independent (at least in our 
numerical precision) which is a behavior completely dif- 
ferent from the two dimensional spin glass (paramagnetic 
phase) . For a paramagnetic phase and very long waiting 
time (i.e. all the points lie in the 1 — C straight line, see 
for example our figure ^ we have equilibrated the sys- 
tem) the PaT scaling plot should consist in points over 
the linear part (quasi-equilibrium regime) , and no one in 
the power law part (aging regime). 

We finally remark that this scaling in addition with the 
analysis of the OFDR (see above) provides us with a pic- 
ture that could be explained assuming a low temperature 
phase with RSB. 

IV. DISCUSSION AND CONCLUSIONS 

We have studied how the fluctuation dissipation rela- 
tions work off equilibrium in the three dimensional Ising 
spin glass with a magnetic field. 

We have shown numerical data that has been obtained 
simulating very large lattices {L = 20, 30 and 60) and for 
extremely large times for the three dimensional Ising spin 



glass. In order to achieve these lattice sizes and time we 
have used a dedicated machine (SUE). 

We have identified with this tool a paramagnetic phase 
in the high temperature region (as expected) and a phase 
where we have found strong violations of fluctuation- 
dissipation. We can describe very well (within our statis- 
tical precision) these violations assuming a RSB scenario, 
yet we can not exclude completely a droplet scenario. 

Moreover we have shown the crossover, both moving 
in temperature as well as moving in magnetic field, be- 
tween a spin glass behavior and a paramagnetic one. We 
have a picture of the phase diagram composed by three 
points (using the notation {T,h)): (1.138,0), (-- 1.1,0.2) 
and (0.714, ~ 0.6) (with the symbol ^ we denote that 
the figure that follows is only indicative). In addition 
we know that for the Gaussian model there is a critical 
point at (0, 0.6), but this critical magnetic field should be 
modified to take into account that we are using binary 
couplings. 

Finally we have checked that the overlap probability 
distribution of our model, P{q), (obtained via the static- 
dynamics link: X{C) ^il)), which is not trivial (at 
least within our numerical precision, see section II.), sat- 
isfies the PaT Ansatz. 
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